Suppressing the Rayleigh- Taylor instability with a rotating magnetic field 
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The Rayleigh- Taylor instability of a magnetic fluid superimposed on a non-magnetic liquid of 
lower density may be suppressed with the help of a spatially homogeneous magnetic field rotating in 
the plane of the undisturbed interface. Starting from the complete set of Navier-Stokes equations for 
both liquids a Floquet analysis is performed which consistently takes into account the viscosities of 
the fluids. Using experimentally relevant values of the parameters we suggest to use this stabilization 
mechanism to provide controlled initial conditions for an experimental investigation of the Rayleigh- 
Taylor instability. 

PACS numbers: 75.50.Mm,47.20.Ma 
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I. INTRODUCTION 

The Rayleigh-Taylor instability Q, 0, is a classi- 
cal hydrodynamic instability Q with relevance in such 
diverse fields as plasma physics, astrophysics, meteorol- 
ogy, geophysics, inertial confinement fusion and granular 
media, for a review see, e.g., 0. Generically this insta- 
bility develops if a layer of liquid is superimposed to an 
immiscible and less dense liquid such that the potential 
energy of the system can be reduced by interchanging 
the liquids. Consequently, the initially plane interface 
between the liquids becomes unstable and the charac- 
teristic dimples and spikes develop resulting finally in a 
stable layering with the lighter fluid on top of the heavier 
one. 

A quantitative experimental investigation of the 
Rayleigh-Taylor instability requires reliable control of the 
initial condition. Standard procedures like suddenly re- 
moved partitions between the fluids 0, or quickly 
turning the experimental cell upside down |3] , clearly 
produce unpredictable initial perturbations. It is much 
more convenient to use some additional mechanism which 
first stabilizes the unstable layering of the liquids and may 
later be switched off instantaneously. It is well known 
that the Rayleigh-Taylor instability may be supp ressed, 
e.g., by vertical oscillation of the system 0], (llj and by 
appropriate temperature gradients |l2j| . The first mech- 
anism is likely to induce uncontrolled initial surface de- 
flections when stopped, in the second one it is difficult to 
abruptly switch off the stabilization. 

In the present paper we investigate the possibility to 
stabilize a potentially Rayleigh-Taylor unstable system 
involving a magnetic fluid by external magnetic fields. 
We will show that for experimentally relevant parame- 
ter values moderate fields strengths which can easily be 
switched on and off are sufficient to achieve the desired 
stabilization. 

Magnetic fluids are suspensions of ferromagnetic nano- 
particles in carrier liquids with the hydrodynamic prop- 
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erties of Newtonian fluids and the magnetic properties 
of super-paramagnets 0] , 0] . It is well known that a 
magnetic field parallel to the plane interface between a 
ferrofluid and a non-magnetic fluid suppresses interface 
deflections with wave vector in the direction of the field 
[l3j . This may be used to stabilize the Rayleigh-Taylor 
instability in two-dimensional situations where the inter- 
face is line, as e. g. in a Hele-Shaw cell. In the more nat- 
ural three-dimensional setting we are interested in here a 
static magnetic field parallel to the undisturbed interface 
is not sufficient to stabilize the flat interface since pertur- 
bation with wave vectors perpendicular to the magnetic 
field will still grow as in the absence of the field. We 
therefore propose to use a spatially homogeneous mag- 
netic field rotating in the plane of the undisturbed inter- 
face and determine appropriate values of the field am- 
plitude and rotation frequency An alternative possibil- 
ity is to use a static inhomogeneous magnetic field with 
the magnetic force counterbalancing gravity |17|. This 
method was used in 0] to investigate the 2-d Rayleigh- 
Taylor instability in a Hele-Shaw cell. 

Before embarking upon the detailed analysis we would 
like to mention three characteristic features of our 
method. Firstly, we will show that a rotating magnetic 
field is unable to suppress all possible unstable modes of 
the system. In fact it can only stabilize surface deflec- 
tions with wavenumber modulus larger than some thresh- 
old value. Perturbations with very long wavelength arc, 
however, not a serious problem in real experiments be- 
cause these are suppressed automatically by the finite ge- 
ometry of the sample. Secondly, it is well-known that in 
analogy with the Faraday instability an oscillating mag- 
netic field will induce new instabilities at wave numbers 
which without field were stable [2(j, |2l]. In our analy- 
sis we keep track of these unstable modes and determine 
the magnetic field strength such that no new instabilities 
may occur. For suppressing these new modes viscous 
losses in the liquids will be decisive which is the reason 
why the viscosities of the two liquids will be consistently 
taken into account in the analysis. Finally, due to the 
dispersed magnetic grains ferrofluids have usually com- 
paratively high densities. We therefore specialize to the 
case in which the upper, heavier layer is formed by the 
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magnetic fluid. This should be the typical situation in 
experiments. Nevertheless a similar analysis with analo- 
gous results is possible for the reverse situation with the 
ferrofluid at the bottom of the system superimposed by 
an even denser non-magnetic liquid. 

The paper is organized as follows. In section II we col- 
lect the basic equations and boundary conditions. In Sec- 
tion III we linearize these equations around the reference 
state of a plane interface between the liquids. Section IV 
contains the Floquet theory to determine the boundaries 
separating stable from unstable regions in the parameter 
plane. After shortly discussing two approximate treat- 
ments of the fluid viscosities in section V we present the 
results of our analysis in section VI. Finally section VII 
contains some discussion. 



II. BASIC EQUATIONS 



We consider a ferrofluid with density p^ superimposed 



on a non-magnetic fluid of lower density p ( ' < p ( ' , see 
Fig. ^ Both layers are assumed to be infinite in hori- 
zontal as well as in vertical direction. The densities and 
the respective viscosities rj^ and rf 2 ^ are taken to be 
constant. The liquids are immiscible and the interface 
between them is parametrized by z = £(x,y,t). We will 
study the stability of a flat interface which we take as 
the rr-y-plane of our coordinate system, the undisturbed 
interface is hence given by £(x,y,t) = 0. In the absence 
of a magnetic field this situation is unstable due to the 
Rayleigh-Taylor instability HHI3. 

In the presence of an external magnetic field H the 
magnetic fluid builds up a magnetization M which is as- 
sumed to be a linear function of the field, M = %H, 
where \ denotes the susceptibility related to the relative 
permeability by p r = 1 + \. Both liquids are subject 
to the homogeneous gravitational field acting in negative 
z-direction and to the interface tension a acting at their 
common interface. The magnetic fluid is additionally in- 
fluenced by the magnetic force density (MV)H result- 
ing from the externally imposed spatially homogeneous 
magnetic field Ho = Ho(cos(ilt) , sin(fii), 0) rotating with 
constant angular frequency f£ in the x-y-planc. 

The time evolution of the system is governed by the 
following set of equations. The incomprcssibility of both 
liquids gives rise to the continuity equations for the ve- 
locity fields v") 



V • v (i) = , 



(1) 



with j — 1,2 where here and in the following the lower 
(i.e. non-magnetic) fluid parameters are denoted with 
superscript (1) and the higher (magnetic) one with su- 
perscript (2). The hydrodynamic equations of motion are 
the Navier-Stokes equations 



p u) (d t + • v) = V • + p {i \ 



(2) 



z ' 














FIG. 1: Sketch of the system under consideration. A fer- 
rofluid of relative permeability p r , density p^ 1 ' and viscosity 
rf 2 ' is superimposed on a non-magnetic fluid with density 
p^ 1 < p^ and viscosity r) '. The normal vector n on the 
interface £(x, y, t) points into the non-magnetic liquid. The 
vector g denotes the gravitational acceleration, a is the inter- 
face tension. The system is subjected to a horizontal rotating 
magnetic field Ho = ffo(cos(Sli), sin(f2i), 0) with angular fre- 
quency SI. 



with g = (0, 0, — g) denoting the acceleration due to grav- 
ity and the stress tensors given by 



T 



(3) 



Here pv> denotes the pressure in each liquid, and B^' = 
Pq(H^ + M^) is the respective magnetic induction. 
Note that also the stress tensor for the non-magnetic liq- 
uid contains contributions from the magnetic field which, 
however, are divergence free and therefore do not give rise 
to a force density in the lower fluid. 

For values of fl relevant to the present investigation 
radiative effects are negligible and the magnetic field has 
to obey the magneto-static Maxwell equations 



v • = 

V x H w = 



(4) 



In view of the second equation it is convenient to intro- 
duce scalar magnetic potentials <f>^ according to = 
— V<&^. These potentials then fulfill the Laplace equa- 
tions 



V 2 $ (J) = 



(5) 



The above equation have to be supplemented by appro- 
priate boundary conditions. Far from the interface the 
velocities must remain bounded, 



lim 

z — >±oc 



< OO 



(6) 



and the magnetic field must be equal to the externally 
imposed field, 



lim -V$ 0) = H G 

z— »±oo 



(7) 



3 



To formulate the boundary conditions at the interface 
z = ((x, y, t) we define the normal vector n by 



V(z - C(x,z,t)) 



(8) 



\V(z-t(x,z,t))\ ■ 
The normal component of the stress tensor has to fulfill 



n e AT e 



-oK 



(9) 



where K = V • n is the local curvature of the interface 
and the symbol AA = A' 2 ' — A^ 1 ) denotes here and in 
the following the difference in the value of the respective 
quantity slightly above and slightly below the interface. 
The tangential components of the stress tensor have to 
be continuous. 



tf AT em n m = 



(10) 



for all vectors t perpendicular to n. The motion of the 
interface is related to the velocity fields in the liquids by 
the kinematic condition 



(11) 



Finally, at the interface the normal component of B and 
the tangential component of H have to be continuous 
which gives rise to the following boundary conditions for 
the magnetic potentials at the interface 



(12) 



$(!) = $( 2 ) 

III. LINEAR STABILITY ANALYSIS 



The main purpose of the present work is to investi- 
gate whether the Rayleigh- Taylor instability can be sup- 
pressed with the help of a rotating magnetic field. We 
will hence study the linear stability of the reference state 
with a flat interface, £(2;, y, t) = 0, in dependence on the 
magnetic field strength Hq and the angular frequency fl. 
The reference solution of the basic equations is given by 



V« = 0, j#> = -pWgz, 

$ = -H (cos(fii)a: + sm{Vlt)y) 



(13) 



To investigate its stability we introduce as usual small 
perturbations 



(14) 



and linearize the basic equations in these perturbations 
as well as in the interface deflection ((x,y,t). We will 
denote the components of the perturbed velocity vectors 



by Sv^' = («w,o®,ti)»'). It is convenient to introduce 
dimcnsionlcss quantities according to 



k c 



t 

p(j) 
H 2 



g*Ap 



1/4 



\J Apga p^' 
( 2 x + 2 



\Po X 2 
Ap p& 



yjApga H 



(15) 



V 



Apa 1 



3\ 1/4 



,0') 



where k c ^H=a = \J Ap g /a is the critical wave number of 
the Rayleigh- Taylor instability for Ho = 0. 

The linearized set of basic equations (|1I2I5[) then reads 



V • <5v w = 

_ 2 x±2 v (h£>.v^: 

X v 

S7 2 d> {l) = . 



(16) 
(17) 
(18) 

(19) 



In order to eliminated the pressure it is convenient to 
consider the z-componcnt of the curl curl of eq. I|17|) which 
is of the form 



(d t -^W 2 ) v 2 w^ = o 



(20) 



where we have introduced the kinematic viscosities v^> 
From the boundary conditions JBJ) and Q we find 



and 



lim io W < oo 

z — >±oo 



lim d z <j) [j) = 

z — >±oo 



(21) 



(22) 



The boundary conditions at the interface simplify under 
linearization. Generally we may replace the interface po- 
sition z = y, t) by z — to linear order in (,{x, y, t). 
Therefore the symbol A has now the more specific mean- 
ing A A = limjijo A(z) — lim z |o -M 2 )- From (|11|) we then 
get 



dtC = w 



U) 



2=0 



implying 



Aii) = 



(23) 



(24) 
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Moreover, the continuity of the flow field v together with 
(|TT))) gives rise to 



From (U3 we find 



A 



A d r w = 



V (Vi - dl) 



(25) 



= 



(26) 



Finally, linearization of © together with H17|) yields 



A 



(M-?7(3Vi + 9 z 2 ))^ 



C + 2- H - V<; 



A 



(27) 



where the horizontal Laplace operator is defined by = 

d 2 + d 2 . 

The magnetic boundary conditions i|12[l acquire the 
form 

A0 = 

M ( cos(nt)d x ( + sin(fii)5 y c) + d z (/^</> (2) - </> (1) ) = . 

(28) 

To find a solution of the set of linearized equations 
(|T$)l . lj2l7|l together with their boundary conditions we 
may exploit their translational invariance and have to 
keep in mind their explicit time dependence induced by 
the second boundary condition ij^S)! for the magnetic field 
problem. An appropriate ansatz is therefore given by 



w^(x,y,z,t) 



and 



(1) (x,y,Z,t) = 4>(t) e ^x+k v y)+kz 

4> {2) {x,y lZl t) = $(t) e l ^ x+k y y ^- kz 



(29) 



(30) 
(31) 



With the abbreviation k 2 
the form 



k 2 eq. (|20() acquires 



(d t - i/CJ") (d 2 z - k 2 ) ) (d 2 z - k 2 ) (z, t) = 0. (32) 

Moreover the ansatzes (|3U|I and (|31|) already fulfill 1)190. 
(|22|l and the first of the boundary conditions H28|) . The 
second one yields 

kt) = i-X—Ho 1 ^- (cos(nt)k x + sm(nt)ky) , (33) 
X + 2 k 



which gives rise to 
H • S7<t (2) 



z=0 



—X-^H 2 k cos 2 (nt)((t)e l{k * x+k y y) . 



The boundary conditions I|21I23I27|I assume the form 



lim w^' < oo 

z — >±oo 



dtC = w (3) 
Aii = 
A d z w = 



z=0 



A 



V (k 2 +d 2 z ) 



w 







(35) 

(36) 
(37) 
(38) 

(39) 



and, using also l|H4|l. 



A 



pd f - v (d 2 z - 3k 2 )) d z 

-l + H$(l + cos(2ftt))fc + k 2 ) k 2 C . 



(40) 



We now invoke Floquet theory [l8L ITgj j to solve this sys- 
tem of linear differential equations with time periodic 
boundary conditions for the amplitudes w^'(z,t) and 

C(*)- 



IV. FLOQUET THEORY 

In order to analyze the stability of the flat interface we 
employ the following Floquet ansatz for the time depen- 
dence of the interface perturbation amplitude £ and the 
2-components of the velocity w^' : 



2in Qt 



{m™ {j \z,t)} = e^ nt j2{tn,^\z)}< 

n 

(41) 

where a + i(3 is the Floquet exponent. Here a is a real 
number and negative a describes stable situation whereas 
positive a signals an instability of the reference state. 
The imaginary part (3 of the Floquet exponent is cither 
zero or one and distinguishes between harmonic {(3 = 0) 
and subharmonic (/3 = 1) response of the system [Lsj . 
Plugging (jUJ) into (J2U) we find 

(d 2 z - q ^ 2 )(d 2 z -k 2 )wlp(z)=0, (42) 

where 



3£° = • 
Ea. (|4*2*|) has the solution 



a + i(/3 + 2n) 
1M 



n 



(43) 



(34) 



w& ( z ) = A$ e kz + Blp e- kz + D$p e 9 - >z + e" 9 - >z , 

(44) 

where the constants An ...D^ can be determined with 
the help of the boundary conditions I|35I39|I . As a result 
the amplitude of the z-componcnt of the velocity may be 
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expressed in terms of the interface amplitude Q n accord- where 
ing to 



(45) 



B n = 



C 



D, 



'k - <£ } ) (V 2 ) (k + qi 2) ) + (k + 
A v k 3 + (r,UqP + ^qV) fc 2 + („M (fc + <£>) + „W (fc + q {2) )) q£ )2 + 2^q£\£h 

v {2) (fc + gi 2) ) +V {1) + 
( ? / 1 )fc + ^) 9 i 2) )U 2)2 -fc 2 ) 

2fct/ 2) 



,(2) 



,(2) 



(fc - (?/ 2 > (fc + q { n } ) + (k + q£ 

U^q^+^k) (k + qW) 

' 2fct/ 2) 



11^ (k + q^j +r/(D (k + <£ ) 
Finally, using these results in l|40|) we find a relation of the form 

°° ( r 1 

^ |W„Cn-#0 Cn + ^CCn-l+Cn+l) 



fc 3 e (a+i(p+2n)) Ut _ q 



(46) 



(47) 



Since this equation has to hold for all values of t all coefficients in the curly brackets must vanish separately. We 
therefore end up with an infinite homogeneous system of linear equations for the amplitudes £ n in which the off- 
diagonal terms arise due to the time dependence in (j40(l . Nontrivial solutions for the £„ require that the determinant 
of the coefficient matrix vanishes which yields the desired relation between the parameters of the problem. Hq, k 
and a. In the present investigation we are mainly interested in the stability boundaries in the parameter plane. We 
therefore specialize to the case a = and find for the coefficients W n in 147|) 



W n = 



-i(2n + p) fi (pM(kB + qV>D) + p (1 \kA + q^C)j + v ^(k 3 B + q^ 3 D) + v ^(k 3 A + q^ 3 C) 
- 3 (pW(kB + q&D) + 7? (1) (fci + gi 1} C)) + 1 - t 



(48) 



To exploit the solvability condition for a numerical de- 
termination of the stability boundaries we have to trun- 
cate the infinite system of linear equations at some fi- 
nite value n max of n. Comparing the results for different 
values of n max the accuracy of the procedure may be 
estimated. For the results presented in section VI we 
have used 7i max = 19, i.e. we have included 39 terms, 
-19 < n < 19. 



V. SPECIAL CASES 



Before presenting explicit results of our analysis for ex- 
perimentally relevant parameter values it is instructive 
to consider two limiting cases for which alternative ap- 
proaches are available. Let us first discuss the situation 
of ideal liquids, = r/ 2 ) = 0. Using J32J), © and 
(|37J) we may then express w^' in terms of Plugging 
the result into the boundary condition l|4(J|l we obtain 
the following Mathieu equation for the amplitude of the 
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surface deflection £(£): 

c 



Pi + P2 



k + Hlk 2 + k 3 



H^k 2 cos(2fii) = 



(49) 



From the standard stability chart of the Mathieu equa- 
tion we are now able to determine the threshold for 
the amplitude Hq of the external field necessary to sta- 
bilize interface deflections with wavenumber modulus k. 
However, since most ferrofluids are rather viscous this 
theory will not adequately describe the experimental sit- 
uation. 

It is possible to approximately incorporate the influ- 
ence of viscosity by assuming that the dominant contribu- 
tion to viscous damping originates far from the interface 
in the bulk of the fluids where the flow field is identical 
to the one of ideal liquids [U . One may then derive 
a damped Mathieu equation for the amplitudes of the 
interface deflection of the form 



d 2 c + 2 7 a t c + 



-k + H 2 k 2 + k 3 



Pi + P2^ (50) 
+ Elk 2 cos(2fii)) = , 



where the damping constant 7 is given by 



7 : 



2^±^k 2 
Pi + Pi 



(51) 



Since the damped Mathieu equation may be mapped on 
the undamped one |22l | we may again employ the stability 
chart of the Mathieu equation to discuss the stabilization 
of a surface deflection mode with wave vector modulus k 
in the presence of small damping. In the following section 
we compare the results of these approximate estimates 
with those of our complete treatment. 



VI. RESULTS 

In this section we display detailed results of our analy- 
sis for a typical experimental combination of a ferrofluid 
and an immiscible non-magnetic fluid which has been 
used in a related experimental investigation [l5j . The 

1690 kgm~ 3 , 

3 „(1) - 



fluid parameters are as follows: p^ 1 ' 



if' = 0.14 Pas, x 



2.6, = 800 kgnT", ^ 



(2) 

0.005 Pas, and a = 0.012 N m _1 . For the capillary length 
A c = 2i:/k c ,H=o we then obtain A c ~ 7 mm. The dimen- 
sionlcss magnetic field amplitude Hq = 1 corresponds to 
a field of Hq = 3.3kA/m, ft = 1 corresponds to a field 
rotating with frequency / = 14.6 Hz. 

In Figs.[21and[3]we show the regions of instability of the 
flat interface in the k-Ho plane. For Hq = all pertur- 
bations are unstable for which the modulus of the wave 
vector is smaller than 1 (in our dimcnsionlcss units, cf. 




A: 



FIG. 2: Stability chart in the plane spanned by the dimen- 
sionless wavenumber k and the magnetic field amplitude Ho 
for the flat interface between a ferrofluid and a non-magnetic 
fluid in a rotating magnetic field. The dimensionless angular 
frequency of the magnetic field is Q — 0.69 corresponding to 
/ = 10 Hz, the remaining parameters are as given in the main 
text. The full lines separating the white regions of stability 
from the gray regions of unstable combinations derive from 
a numerical solution of 1471 . For comparison also the results 
of the inviscid theory building on (I49H and of the approxi- 
mate treatment of viscosity related to 150H are included as 
long and short dashed lines respectively. The threshold val- 
ues of the magnetic field and wave number are Hth = 2.2 and 
kth = 0.2 corresponding in physical units to H t h — 7.3 kAm -1 
and \ t h ~ 3.7 cm. 




h. 



k 



FIG. 3: Same as Fig. [21 for an angular frequency of SI = 1.37 
(/ = 20 Hz). Due to the higher frequency the influence of 
the viscosity is stronger giving rise to the threshold values 



Hth =3.1 and kth 
and A t h — 7.4 cm. 



0.1 corresponding to H t h — 10.3 kAm 



(|15fl h which is the well-known trademark of the Rayleigh- 
Taylor instability. Increasing H from zero the interval of 
unstable wave numbers shrinks and hence more and more 
long-wave perturbations may be stabilized. However, if 
Hq gets larger than a threshold value H t h the parametric 
excitation due to the time-dependent magnetic field gives 
rise to new instabilities at higher wave-numbers. Since 
these additional unstable modes are clearly unwanted 
Hq must remain below this threshold value H t h- Cor- 
respondingly there is a threshold kth f° r the wavenumber 
modulus such that perturbations with k < k t h cannot be 
stabilized with the help of the magnetic field. As we will 
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detail in section VII these modes have to be stabilized by 
lateral boundary conditions. We note that with decreas- 
ing n the tongues of instability move closer together and 
come nearer to the fc-axis implying H t h — * and k t h — > 1 
for Q -> 0. 

It is clearly seen from the figures that the stability 
regions are strongly influenced by the viscosity of the 
liquids. In the inviscid theory the tongues of instability 
all reach the fc-axis implying that any rotating magnetic 
field would induce new unstable modes at values of k 
that were stable in the absence of the field. Therefore a 
complete suppression of the Rayleigh- Taylor instability 
would be impossible. It is also apparent that for realistic 
parameter combinations the phenomenological inclusion 
of viscosity in the theoretical description as discussed in 
the previous section may give results which significantly 
differ from the complete theory. This is similar to the 
analysis of the Faraday instability performed in |19| . 

Figs. 01 and [3] display the dependence of the threshold 
values H t h and kth on the angular frequency O of the 
field. Clearly H th increases and k th decreases with in- 
creasing Q as exemplified also by a comparison between 
Fig. |21 and [3J For the parameters considered an increase 
in f2 beyond Q = 2 docs not significantly reduce kth any 
more. 




FIG. 4: Threshold value Hth of the magnetic field amplitude 
Ho as function of the angular frequency of the field for the 
parameters given in the main text. The displayed interval 
f2 = 1...7 corresponds in physical units to / ~ 15 Hz. ..100 Hz. 
The inset shows a blow-up of the steep increase of Hth for 
small values of Q. 

Finally, Fig. H3 combines Figs. 0] and [5] and shows the 
relation between the two threshold values H t h and kth- 

VII. DISCUSSION 

In the present paper we have investigated the possi- 
bility to stabilize a layering of a ferrofluid and a non- 
magnetic fluid which were potentially unstable due to 
the Rayleigh- Taylor instability by a spatially homoge- 
neous magnetic field rotating in the plane of the undis- 
turbed interface. Special emphasis was put on an exact 
treatment of the influence of the viscosities by starting 
from the complete set of Navicr-Stokcs equations for both 
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FIG. 5: Similarly to Fig. [I] the threshold value k t h of the 
wavenumber modulus k is shown as function of the angular 
frequency Q of the field. 



kth 




6 Ht h 7 



FIG. 6: Threshold value k t h of the wave number k versus 
threshold value H t h of the magnetic field amplitude Ho, again 
for the special set of parameters given in the main text. 



liquids. Our results show that this approach is for exper- 
imentally relevant parameter values superior to both the 
inviscid theory and to a standard phenomenological pro- 
cedure to include viscous effects using the inviscid flow 
field. 

The trademark of the Rayleigh- Taylor instability is a 
band of unstable wave numbers extending from k = 
up to a threshold value kth which in the absence of 
magneti c effect s is given by the capillary wavelength 
k c = \J Apg/a. The main result of the present inves- 
tigation is that k t h may be reduced roughly by a factor 
of ten with the help of a rotating magnetic field of ex- 
perimentally easily accessible amplitude and frequency. 
As expected the stabilization works best for ferrofluids 
with high susceptibility \ which, however, have also high 
densities increasing in turn k c . 

In order to provide a clean initial condition for an ex- 
perimental study of the onset of the Rayleigh- Taylor in- 
stability one has also to stabilize the modes with k < kth- 
One way to accomplish this suppression may be to use 
the boundary condition of a finite geometry, i.e. by pin- 
ning the contact line between the liquids at the boundary 
of the sample. In this way all long wave-number modes 
with k < kh c are stabilized. Here k(, c is determined by 
the linear extension L of the sample and roughly given 
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by kbc — tt/L. Modes with k > k c ,H=o arc suppressed 
by surface tension. If one is able to temporarily stabilize 
the remaining modes by the rotating magnetic field, i.e. 
if one is able to realize k t h < kb c the flat interface is sta- 
ble. Switching off the magnetic field at a given time all 
modes with kbc < k < 1 will become unstable. Since it is 
easily possible to realize values of k t h significantly smaller 
than k max , the wavenumber with largest growth rate in 
the absence of the field, the ensuing Rayleigh- Taylor in- 
stability should closely resemble the case without lateral 
boundary conditions. To be precise it should be empha- 
sized that our theoretical analysis is for infinite layers 
only and does not take into account the influence of lat- 
eral boundary conditions. However, the relevant values 
of Hth and kth will only marginally be modified. 

For an order of magnitude estimate let us consider a 
cylindrical vessel of diameter d = 5 cm. Pinning the con- 
tact line at the boundary the instability of modes with 



dimcnsionlcss wavenumber k < kb c — 0.11 will be sup- 
pressed. On the other hand a rotating magnetic field with 
amplitude Hq ~ lOkArn -1 and frequency / = 20 Hz re- 
alizes kth — 0.1 (cf. Fig- EI) • Switching off the magnetic 
field all modes with kb c < k < 1 will become unstable. 
For the above example this includes the first 8 cylindrical 
modes which should allow a rather accurate study of the 
Rayleigh- Taylor instability. We hope that our theoreti- 
cal study may stimulate experimental work along these 
lines. 
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